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Abstract 


Much  work  is  currently  being  carried  out  in  the  field  of  mathematical 
modelling  of  the  physical  and  chemical  processes  taking  place  within 
combustion  equipment  to  enable  better  combustion  chamber  design.  One 
aspect  of  this  is  the  numerical  prediction  of  burning  droplet  trajectories 
in  a  gas  flow  field  and  in  this  paper  a  general  three-dimensional  particle/ 
droplet  trajectory  algorithm  is  presented  and  applied  to  fuel  sprays  in 
gas  turbine  combustors.  Droplet  evaporation  and  initial  heating  up 
effects  are  included  and  trajectories  are  calculated  within  a  three- 
dimensional  computer  predicted  gas  flow  field.  By  injecting  many 
droplets  with  initial  conditions  representing  a  particular  fuel  nozzle 
a  complete  fuel  spray  can  be  constructed,  and  this  is  presented 
graphically.  The  results  show  that  in  some  cases  the  three-dimensional 
nature  of  the  gas  flow  and  the  behaviour  of  large  droplets  may  be  expected 
to  have  a  significant  effect  on  combustor  performance. 


-  1  - 
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Evaporation  constant 
Drag  coefficient 

Specific  heat  capacity  at  constant  pressure 
Droplet  or  particle  diameter 
Rosin-Rammler  mean  diameter 

Finite  difference  grid  subscripts  in  axial,  radial  and 
tangential  directions. 

Latent  heat  of  evaporation 

Mass 

Mass  flowrate 

Rosin-Rammler  spread  parameter 

Nusselt  number 

Prandtl  number 

Reynolds  number 

Time 

Temperature 
Velocity  vector 

velocity  components  in  axial,  radial  and  tangential  directions 
time  derivatives  of  u,  v  and  w 

distance  in  axial,  radial  and  angular  directions, 
time  derivatives  of  x,  y  and  9 
Densi ty 
Viscosity 

Thermal  conductivity 


Subscripts : 


o 

CO 

g 

l 

p 


Initial 
Gas  flow 
Gas 

Liquid 

particle  or  droplet 
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Introduction 


The  availability  of  large  general  purpose  digital  computers  has 
enabled  numerical  techniques  to  be  applied  to  many  complex  engineering 
problems.  In  the  field  of  Combustion  Research  mathematical  modelling 
can  be  used  to  predict  the  behaviour  of  all  kinds  of  combustion  processes 
and  such  numerical  procedures  are  forming  the  basis  of  a  more  systematic 
approach  to  combustion  equipment  design.  The  aim  of  current  work  on 
Gas  Turbine  Combustor  modelling  is  to  provide  a  design  tool  to  assist 
in  what  is  at  present  a  largely  empirical  process  of  combustor  design. 

This  is  a  field  of  increasing  importance  as  a  result  of  present  interest 
in  fuel  economy  and  pollution  control. 

Algorithms  have  been  produced  capable  of  modelling  the  chemical  and 
physical  processes  occuring  in  a  gas  turbine  combustor  can  and  these 
are  reviewed  by  Swithenbank  et  al  (Ref.  1).  These  models  do  not,  however, 
include  the  details  of  two-phase  flow  and  droplet  combustion  arising 
from  the  injection  of  a  liquid  fuel  spray.  The  fate  of  burning  droplets 
in  a  gas  flow  has  so  far  only  been  considered  in  two  dimensions  (Ref.  2). 
Other  related  work  includes  that  of  Law  (Ref.  3)  who  analyses  the  motion  of 
an  evaporating  droplet  in  a  constant  cross  flow,  and  Chiu  et  al  (Ref.  4) 
who  present  a  model  for  the  collective  behaviour  of  burning  droplets  in 
liquid  fuel  sprays.  Mellor  (Ref.  5)  has  represented  the  dominant  processes 
of  spray  combustion  by  means  of  time  scales  for  evaporation,  mixing 
and  chemical  reaction.  Dring  and  Suo  (Ref.  6)  have  obtained  particle 
trajectories  in  swirling  flows  and  subsequently  in  turbine  cascades  (Ref.  7). 

At  present  the  detailed  behaviour  of  fuel  droplets  during  spray 
combustion  is  largely  unknown  and  experimental  data  are  difficult, 
if  not  impossible  to  obtain.  These  phenomena  affect  all  aspects  of  the 
combustion  process  and  a  full  understanding  of  spray  combustion  is 
required  before  rigorous  design  techniques  can  be  employed. 

In  this  study  a  general  algorithm  is  produced  to  calculate  the 
trajectories  of  particles  or  evaporating/buming  droplets  in  a  three- 
dimensional  gas  flow  field  represented  by  a  finite  difference  grid,  and 
this  is  applied  to  fuel  sprays  in  a  gas  turbine  combustor.  The  results 
predict  in  detail  the  behaviour  of  individual  burning  droplets  and  the 
distribution  of  vaporised  (unburnt)  fuel  in  the  combustor  can. 


Theory 


Because  of  the  cylindrical  nature  of  combustion  chambers  and  the 
gas  flows  involved,  a  cylindrical  polar  coordinate  system  is  used  to 
describe  the  location  of  the  particle  or  droplet.  The  equations  of 
motion  of  a  particle,  neglecting  all  forces  except  drag,  in  component 
form  are: 
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where  u  ,  v  and  w  are  the  absolute  particle  velocities  in  the 
P  P  P 

radial  and  tangential  directions  respectively,  and  u^,  v^and 
the  corresponding  gas  velocities. 

For  the  cylindrical  polar  coordinate  system,  the  equations 


axial , 
are 

of  trajectory 


are : 


in  the  axial,  radial  and  angular  directions. 

In  the  above  equations  the  relative  Reynolds  number,  Rg  is  defined 

as : 
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where  Up  and  u^  are  the  velocity  vectors  of  the  particle  and  gas  stream 

respectively. 

The  drag  coefficient,  C^,  is  given  by  the  following  equations 
(Ref.  8): 
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CD  =  27  Re~°‘84  0<Re<80 


C  =  0.271  Re°'217  80<Re <104 
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Equations  (1)  to  (10)  describe  the  trajectory  of  a  particle 
in  a  gas  stream.  For  an  evaporating  droplet,  is  not  constant, 
and  the  rate  of  change  of  droplet  diameter  with  time  is  required. 
For  forced  convection  this  is  given  by: 


dD^  r 
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dt  ^  (1  +  0,23  ^  (11) 

P 

where  C  is  the  evaporation  constant.  This  is  dependant  on  the  properties 
b 

of  the  fuel  as  well  as  the  surrounding  gases,  and  a  widely  used  expression 

for  C  derived  from  the  quasi-steady  analysis  of  droplet  combustion  is 
b 


that  given  by  Wise  and  Agoston  (Ref.  9) 
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To  allow  for  droplets  "heating  up"  after  entry  into  the  hot  gas 
stream  it  has  been  assumed  that  droplets  do  not  begin  to  evaporate 
until  they  reach  their  boiling  point.  The  rate  of  change  of  temperature 
with  time  is  given  by: 
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where  T  is  the  droplet  temperature. 

P 

If  it  is  assumed  that  the  Reynolds  number  varies  slowly  with  respect 
to  time,  it  is  possible  to  integrate  equations  (11)  and  (13)  for  small 
increments  of  time  so  that  equation  (11)  becomes: 
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giving  the  droplet  diameter  at  time  t,  where  D^is  the  droplet 

diameter  at  time  t  .  Equation  (13)  is  integrated  to  give: 
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where  Tp0  is  the  drop  temperature  at  time  tQ. 

Equation  (13).  is  suppressed  until  the  droplet  temperature, 

reaches  the  boiling  point  of  the  liquid,  after  which  T  remains  constant. 

P 

The  Nusselt  number,  Nu  is  given  in  Refs.  10  and  11  as: 
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where  the  Prandtl  number,  Pr  is  given  by: 

Pr  =  CP8  US  (17) 

Xg 

Equations  (1)  to  (17)  above  describe  the  motion  of  an  individual 
evaporating  droplet  in  three  dimensions  in  a  gas  stream. 

Numerical  integration  can  now  be  used  to  calculate  the  droplet 
trajectory  in  the  given  gas  flow  field. 

Application 

For  the  purposes  of  this  investigation  a  computer  predicted  flow 
field  in  a  gas  turbine  combustor  was  used.  This  has  been  calculated 
for  a  Lycoming  combustor  using  a  numerical  prediction  algorithm  to 
model  the  physical  and  chemical  processes  taking  place  (Ref.  12). 

The  geometry  of  the  Lycoming  combustor  can,  which  was  designed 
as  research  combustor,  is  shown  in  Fijg.l.  Air  enters  via  a  swirler  which 
produces  a  finite  swirl  velocity.  The  fuel  spray  is  introduced  at  a 
point  near  the  centre  line  in  front  of  a  baffle  and  is  of  hollow  cone 
type.  Primary,  secondary  and  dilution  air  streams  are  introduced  by  three 
sets  of  injection  holes,  each  consisting  of  six  equally  spaced  circular 
orifices  around  the  periphery  of  the  combustor  can.  This  is  responsible  for 
the  three-dimensional  nature  of  the  problem,  and  the  cyclic  nature  of 
the  flow  enables  the  combustor  can  to  be  represented  by  a  single  60° 
sector,  which  is  divided  into  a  27  x  18  x  7  element  grid  in  the  x,  y 
and  8  directions  respectively.  The  grid  is  shown  diagramatical ly  in  Fig.  2. 

Each  grid  node  is  encompassed  by  a  computational  cell  which  is  used 
to  store  temperature  and  density  data.  Gas  velocity  data,  however,  is 
defined  for  staggered  control  volumes  as  shown  in  Fig.  3.  For  droplet 
trajectory  calculations  these  data  are  redefined  to  refer  to  the  same 
nominal  control  volumes  as  used  to  store  temperature  and  density  data, 
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by  averaging  adjacent  velocity  values.  This  considerably  simplifies 
the  subsequent  trajectory  calculations. 


The  two-dimensional  model  of  Reference  2  assumes  that  the  gas 
flow  field  is  essentially  two-dimensional  in  the  region  of  the  fuel 
spray.  To  investigate  the  validity  of  this  assumption,  a  statistical 
analysis  of  axial  velocity  data  was  made.  For  each  value  of  I  and 
each  of  three  J  values  the  set  of  data  for  all  K  values  was  averaged 
and  sample  standard  deviation  calculated.  The  result  of  this  analysis  is 
given  in  Fig.  4  and  shows  that  the  flow  is  essentially  two-dimensional 
until  a  short  distance  upstream  of  the  first  (primary)  set  of  air 
injection  holes.  Beyond  this  point  the  flow  becomes  predominantly 
three-dimensional,  particularly  in  the  vicinity  of  the  air  injection  holes. 
It  can  also  be  seen  that  flow  is  axisymetric  near  to  the  axis  (J-3). 

A  vector  plot  of  the  velocity  field  is  given  in  Fig.  5.  It  can 
be  clearly  seen  that  there  is  a  recirculation  zone  in  front  of  the  baffle, 
and  that  the  flow  in  the  region  of  the  air  injection  holes  is  dominated 
by  the  inward  flow  of  air.  The  flow  cross-section  at  the  point  of  primary 
air  inlet  shows  how  the  three-dimensional  nature  of  the  flow  field  is 
brought  about  by  the  air  injection. 

Once  the  u-,  v-,  w-velocity,  temperature  and  density  data  is  stored 
as  a  27  x  18  x  7  element  array,  the  fundamental  equations  of  motion  and 
trajectory  (equations  1  to  6)  can  be  solved  numerically  by  assuming  that 
the  gas  velocity,  temperature  and  density  remain  constant  throughout 
each  computational  cell.  These  six  simultaneous  differential  equations 
are  solved  by  using  the  fourth  order  Runge-Kutta  method  for  a  short  time 
step,  the  length  of  which  is  adjusted  automatically  to  give  approximately 
the  same  number  of  time  steps  while  the  droplet  is  in  each  cell. 

After  solving  the  equations  for  one  time  step  the  new  drag  coefficient 
and  Reynolds  number  are  calculated  and  one  of  equations  (14)  or  (15)  is 
used  depending  on  whether  the  droplet  has  reached  its  boiling  point. 
Computation  for  the  next  time  increment  then  proceeds  unless  the  droplet 
has  passed  out  of  the  computational  cell.  In  this  case  the  new  cell  that 
the  droplet  occupies  is  determined  and  the  iteration  is  resumed.  At 
this  stage  a  check  is  also  made  for  droplets  passing  into  an  adjacent 
sector  (in  which  case  they  are  made  to  re-enter  the  original  sector 
at  the  symetrically  corresponding  position),  and  for  other  contingencies 
which  may  terminate  the  iteration  procedure,  e.g.  droplets  vanishing  by 
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evaporation,  colliding  with  the  walls  of  the  combustor  can,  etc. 


The  amount  of  fuel  evaporated  per  unit  time  in  each  cell  is  calculated 
for  each  droplet  as  i t  leaves  the  cell  as: 
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evap 


3  3 

D  -D. 
out  in 


"initial 


f  ue  1 
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where  D.  is  the  droplet  diameter  on  entering  the  cell,  D  the 
in  °  out 

exit  diameter  and  D.  .  .  _  the  droplet  diameter  at  the  start  of  its 
initial  r 

trajectory.  The  total  mass  flowrate  of  fuel  represented  by  the  droplet  is 
mf ue 1 ' 

If  the  droplet  hits  a  wall  the  remaining  fuel  contained  in  the 
droplet  is  deposited  into  the  cell  adjacent  to  the  wall,  representing 
droplets  adhering  to  the  wall  of  the  combustion  chamber.  It  would  be 
possible  here  to  extend  the  model  to  include  droplets  rebounding  off  the 
wall,  shattering  etc. 


To  extend  the  model  to  predict  the  behaviour  of  a  whole  spray, 
a  statistical  size  distribution  is  used  to  represent  the  spray  as  a 
number  of  size  ranges,  each  represented  in  turn  by  a  single  droplet 
diameter,  the  mean  of  the  size  range.  Inherent  in  this  approach  is  the 
assumption  that  an  individual  droplet  trajectory  is  unaffected  by  the 
presence  of  other  droplets  in  the  droplet  cloud.  This  is  valid  if  the 
rate  of  reaction  (i.e.  combustion)  of  fuel  vapour  is  sufficiently 
high  to  prevent  a  significant  partial  pressure  of  unburnt  fuel  vapour 
existing  between  droplets. 


It  is  assumed  that  after  sheet  break-up,  the  fuel  spray  leaving 
the  atomiser  obeys  a  Rosin-Rammler  size  distribution.  In  this  model 
the  mass  fraction  of  fuel  having  a  droplet  diameter  greater  than 
D  is  given  as: 


"D 


=  exp 


I)\n 
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where  D  is  the  Rosin-Rammler  mean  diameter  and  n  is  the  spread 
parameter.  For  this  study  a  value  66  microns  was  chosen  for  D 
and  the  spread  parameter  is  taken  as  2.5.  These  values  were  derived 
from  experimental  data  for  a  typical  nozzle  obtained  by  Pan  (Ref.  13). 

The  mass  fraction  of  fuel  in  a  given  size  range  is  found  by  subtraction  : 
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where  the  size  range  lies  between  and  D2 . 

The  size  ranges  are  chosen  to  be  of  equal  length  and  cover 
droplet  diameters  from  zero  to  a  diameter  below  which  99%  of  the 
mass  of  fuel  falls.  A  total  of  16  size  ranges  was  chosen  to 
adequately  model  the  complete  spray  without  using  excessive 
amounts  of  computer  time.  The  sizes  range  from  3 . 8p  to  118', :  in 
steps  of  7.6 g. 

Single  droplets  with  diameters  equal  to  the  mean  of  each  size 
range  are  injected  into  the  gas  flow  field  with  initial  x  and  y 
coordinates  corresponding  to  the  point  of  sheet  break-up  of  the  spray. 

In  addition,  for  each  size  range,  droplets  are  injected  trom  six 
different  angular  locations  differing  by  10°.  Consequently  the 
trajectories  of  a  total  of  96  fuel  droplets  are  calculated. 

The  atomiser  characteristics  are  specified  by  spray  cone  angle, 
velocity,  distance  from  nozzle  to  sheet  break-up,  and  the  Rosin-Rammler 
parameters.  From  this  and  general  physical  data  the  entire  spray  cone 
is  constructed. 

The  amount  of  fuel  evaporated  in  each  cell  is  summed  for  the 
whole  spray.  The  mass  flowrate  of  fuel  represented  by  each  droplet 
trajectory  is: 

m  .  ...  . 

m..  ,  =  range  .  m  .  C*-l ) 

fuel  - total 

o 

where  m  is  defined  in  equation  (20)  .  Since  each  size  range  is 

range 

represented  by  six  droplets  of  different  initial  angular  position,  the 
factor  of  1/6  occurs.  The  total  mass  flowrate  of  fuel  is  m 

total 

Results  and  Discussion 

Droplet  trajectory  computations  for  complete  sprays  are  presented 
graphically  for  sprays  of  nominal  included  cone  angles  of  45°  and 
80°  at  an  injection  velocity  of  20m/s  and  also  for  a  45°  cone  angle 
and  50  m/s  injection  velocity.  The  injection  velocity  for  the  spray 
is  the  resultant  velocity  in  the  x-y  plane  acquired  by  each  droplet 
as  it  breaks  away  from  the  liquid  sheet  at  the  point  of  sheet  break-up. 

No  initial  tangential  velocity  is  given  to  the  drop’ets.  The  distance 
between  the  atomiser  and  the  point  at  which  this  occurs  has  been  taken  as 
5mm.  The  figures  show  side  and  end  elevations  of  the  combustor  can 
and  droplet  trajectories  have  been  projected  orthographical ly  onto  these. 
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The  large  arrows  indicate  the  air  injection  locations.  Droplet  burn¬ 
out  locations  are  encircled,  whilst  trajectories  ending  without  a  circle 
have  hit  the  combustor  wall.  In  the  case  of  the  latter,  calculations 
are  terminated  a  short  distance  before  the  wall. 

Considering  figure  6  it  is  observed  that  for  these  typical  values 
of  cone  angle  and  velocity  (as  might  be  expected  from  a  pressure-jet 
type  of  atomiser),  the  majority  of  droplets  are,  in  fact,  evaporating 
in  the  two-dimensional  region  of  the  combustor  can  (<60p) .  However,  a 
significant  proportion  impinge  on  the  combustor  wall  and  of  particular 
interest  is  the  portion  of  the  spray  deflected  by  the  primary  air  injection 
stream  and  hitting  the  opposite  wall  of  the  can.  Some  of  the  largest 
droplets  (>100p)  are  also  observed  to  be  surviving  into  the  area  near  the 
dilution  air  injection  having  sufferred  slight  deflections  by  both 
primary  and  secondary  air  jets. 

Figure  9  shows  a  histogram  of  the  proportion  of  fuel  evaporated 
along  the  length  of  the  combustion  chamber  for  the  45°,  20  m/s  case. 
Although  the  relatively  small  proportion  of  fuel  surviving  into  the 
three-dimensional  region  of  the  combustor  can  would  at  first  appear  to 
justify  a  simpler  two-dimensional  model  of  spray  combustion,  it  must  be 
remembered  that  in  a  field  where  very  high  combustion  efficiencies  are 
already  common  small  losses  of  fuel  are  of  increased  importance.  Also, 
in  applications  such  as  aviation,  where  weight  becomes  a  major  factor, 
small  combustion  inefficiencies  become  magnified  in  financial  terms.  In 
terms  of  pollution  emissions,  only  a  very  small  proportion  of  the  fuel 
needs  to  be  involved  to  produce  significant  effects. 

Figures  7  and  8  show  droplet  trajectories  for  two  extreme  cases. 

Fig.  7  shows  a  spray  with  high  velocity.  More  fuel  is  seen  to  be  being 
deposited  on  the  combustor  wall  further  downstream  than  in  the  case  of  the 
20  m/s  spray.  This  would  obviously  constitute  a  bad  choice  of  nozzle 
characteristics.  In  Fig.  8  the  spray  cone  angle  is  clearly  too  large, 
as  a  significant  proportion  (4%)  of  the  fuel  hits  the  chamber  wall. 

Conclusions 

It  is  shown  that  an  analytical  technique  can  be  used  to  predict 
the  details  of  combustion  of  individual  droplets  and  the  fuel  spray 
as  a  whole.  As  a  design  tool  such  an  approach  would  clearly  be  of 
great  help  to  the  combustion  chamber  designer  and  in  the  matching  of 
nozzle  characteristics. 
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The  significance  of  the  behaviour  of  large  droplet  sizes  is 
demonstrated  as  being  a  likely  cause  of  reduced  combustion  efficiency 
and  the  presence  of  pollutants  in  the  exit  gases,  which  may  arise 
from  locally  fuel-rich  areas  and  unburnt  fuel  appearing  in  the  lower 
temperature  regions  of  the  combustor.  An  additional  hazard  is  possible 
corrosion  of  and  damage  to  turbine  blades  due  to  the  presence  of 
pollutants,  soot  and  unburnt  fuel  droplets  in  the  combustor  exit 
gases. 

The  information  obtained  by  means  of  the  computer  prediction  is 
at  present  impossible  to  obtain  by  experimental  measurements.  A  method 
has  yet  to  be  found  whereby  the  computed  results  can  be  verified  experimentally. 
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LYCOMING  COMBUSTOR, 

GRID  ARRANGEMENT  27x18x7  =3402  TOTAL  POINTS 


Fig.  2.  Finite  difference  grid  for  the  Lycoming  combustor  can 


W-AND  V-VELOCITY  CONTROL  VOLUMES  IN 


Fig.  3.  Location  of  u-,  v-  and  w-velocity  control  volumes 


Injection  Velocity :  20  m/s 


Fig.6  Droplet  Trajectories  for  Complete  Spray 
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Injection  velocity:  50  m/s 


Rg.7  Droplet  Trajectories  for  Complete  Spray 


Included  cone  angle:  80° 

Injection  velocity  :  20  m/a 

Fig.  8  Droplet  Trajectories  for  Complete  Spray 


